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Abstract 

We re-examine the predictions of a two-phase model of shadowing in nuclear 
deep-inelastic scattering in light of new NMC data on the x, A and Q 2 depen- 
dence of the ratios of structure functions. The model, which combines vector 
meson dominance at low Q 2 with diffractive Pomeron exchange for large Q 2 , 
agrees with the observed small, but non-zero, slopes in logQ 2 , which indi- 
cate the importance of higher twist shadowing effects in the transition region, 
0.1 < Q 2 < 1 GeV 2 . We note also that the latest E665 data on the deuteron 
to proton ratio suggests the presence of a small amount of shadowing in the 
deuteron. 
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There has recently been a great deal of progress in the experimental study of the A 
dependence of nuclear structure functions in the small-x region where nuclear shadowing 
occurs. At CERN the New Muon Collaboration (NMC) has been able to extract the slopes 
in the four-momentum transfer squared, Q 2 , of the ratios of structure functions of a variety 
of nuclei for x = Q 2 /2Mv down to ~ 0.001 []]-0], where v is the energy transfer to the target 
and M the nucleon mass. Meanwhile, E665 at Fermilab has extended previous measurements 
on the deuteron and several heavier nuclei to x as low asi~ 10~ 5 0-0]. 

These new measurements provide a tremendous challenge to theoretical models of nuclear 
deep-inelastic scattering. Not only do they cover five orders of magnitude in x, but the 
momentum transfer squared also varies from tens of GeV 2 to values as low as 0.05 GeV 2 . 
For the latter one is (to say the least) hard pressed to justify a parton model description. 
Clearly here one is dealing with a transition region between those regimes where hadronic 
and explicit quark treatments are appropriate. While it is unlikely that any single theoretical 
approach will be totally reliable over the entire range of x and Q 2 , it is important that all 
existing models be put to the test with a view to understanding the nature of the transition 
between high and low Q 2 . 

Our aim is to extend earlier calculations of shadowing in deuterium || and heavier nuclei 
|| to compare with these new data. The approach taken is a two-phase model, similar to 
the work of Kwiecinski and Badelek flQ Hl^l - At high virtuality the interaction of the virtual 
photon with the nucleus is most efficiently parametrized in terms of diffractive scattering 
through the double and triple Pomeron. On the other hand, at low virtuality it is most 
natural to apply a vector meson dominance (VMD) model, in which the virtual photon 
interacts with the nucleons via its hadronic structure, namely the p°, uo and <fi mesons. The 
latter contribution vanishes at sufficiently high Q 2 , but in the region of interest here it is in 
fact responsible for the majority of the Q 2 variation. 

For convenience we shall summarize the key ingredients of the calculation, beginning with 
the diffractive component. For this we use Pomeron (IP) exchange between the projectile 
and two or more constituent nucleons to model the interaction of partons from different 
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S^F 2 A (x, Q 2 ) = - \ A dy f IP/A (y) F?{xn>, Q 2 ), (1) 



nucleons within the nucleus. Assuming factorization of the diffractive cross section [ 13| , the 
shadowing correction (per nucleon) to the nuclear structure function F 2 A from iP-exchange is 
written as a convolution of the Pomeron structure function, F 2 /p , with a distribution function 
("flux factor"), fip/A, describing the number density of exchanged Pomerons: 

rA 

Umzn 

where y = x(l + M x /Q 2 ) is the light-cone momentum fraction carried by the Pomeron (Mx 
is the mass of the diffractive hadronic debris), and xjp = x/y is the momentum fraction of 
the Pomeron carried by the struck quark in the Pomeron. The y dependence of fip/A{y) can 
be calculated within Regge theory []8|-|1^,|14]] (for a survey of other definitions of the "flux 
factor" see Ref. ]n|). Within experimental errors, the factorization hypothesis, as well as 
the y dependence of fi P /A{y) (§H3]; are consistent with the results obtained recently by 
the HI [16| and ZEUS |T7[ Collaborations at HERA from observations of large rapidity gap 
events in diffractive ep scattering. These data also confirm previous findings by the UA8 
Collaboration [TS|] that the Pomeron structure function contains both a hard and a soft 
component: F 2 F \xjp,Q 2 ) = F 2 IP( ' ha ' Td \xip J Q 2 ) + F 2 aoft ' (xjp, Q 2 ). The hard component of 
F 2 IP is generated from an explicit qq component of the Pomeron, and has an xjp dependence 
given by xjp(l — xj P ) |T^JT9|1 , in agreement with the recent diffractive data [fT6| P~8|| . The soft 
part, which is driven at small xjp by the triple-Pomeron interaction HlO| , [l9| , has a sea quark- 
like Xip dependence^], with normalization fixed by the triple-Pomeron coupling constant |20[ . 

The dependence of F 2 /p on Q 2 at large Q 2 arises from radiative corrections to the parton 
distributions in the Pomeron |]TT] , pT| , which leads to a weak, logarithmic, Q 2 dependence 
for the shadowing correction S^F^. Alone, the Pomeron contribution to shadowing would 
give a structure function ratio F^/F® that would be almost flat for Q 2 ^ 2 GeV 2 . The low- 
Q 2 extrapolation of the qq component is parametrized by applying a factor Q 2 /(Q 2 + Ql), 
where Qq ~ 0.485 GeV 2 may be interpreted as the inverse size of partons inside the 



1 In Ref. [19| the soft component is associated the qqg Fock component of the 7* wave function. 



3 



virtual photon. For the nucleon sea quark densities relevant for F^ p ^ Boit > we use the recent 
parametrization from Ref. [22|, which includes a low-Q 2 limit consistent with the real photon 
data, in which case the total Pomeron contribution S^ IP ^F A — > as Q 2 — > 0. 

To adequately describe shadowing for small Q 2 requires one to use a higher-twist mech- 
anism, such as vector meson dominance. VMD is empirically based on the observation that 
some aspects of the interaction of photons with hadronic systems resemble purely hadronic 
interactions P3,23|. In terms of QCD this is understood in terms of a coupling of the photon 



to a correlated qq pair of low invariant mass, which may be approximated as a virtual vector 
meson. One can then estimate the amount of shadowing in terms of the multiple scattering 
of the vector meson using Glauber theory. The corresponding correction (per nucleon) to 
the nuclear structure function is: 

5 (v) F A (x C 2 ) - 1 Q2 V M ^ 6avA (2) 

where Sctva is the shadowing correction to the vector meson — nucleus cross section, calcu- 
lated in Ref. ||, fy is the photon — vector meson coupling strength ]23[ (see also and 
My is the vector meson mass. In practice, only the lowest mass vector mesons (V = p°, u, 0) 
are important at low Q 2 . (Inclusion of higher mass vector mesons, including continuum con- 
tributions, leads to so-called generalized vector meson dominance models [pq] .) Usually one 
omits non-diagonal vector meson transitions (VN — > V'N), as these are not expected to be 
large. The vector meson propagators in Eq.(|2|) lead to a strong Q 2 dependence of 5^F A , 
which peaks at Q 2 ~ 1 GeV 2 . For Q 2 — > and fixed x, S^ V 'F A disappears due to the 
vanishing of the total F A . Furthermore, since this is a higher twist effect, shadowing in the 
VMD model dies off quite rapidly between Q 2 ~ 1 and 10 GeV 2 , so that for Q 2 > 10 GeV 2 
it is almost negligible — leaving only the diffractive term, S^ Ips> F A . (Note that at fixed u, 
for decreasing Q 2 the ratio F^/F^ approaches the photoproduction limit.) 

While the asymptotic Q 2 behavior of nuclear shadowing seems clear, there is still con- 
siderable interest in the transition region where the high- and low-Q 2 descriptions merge. 
In practice, this occurs for Q 2 between about 0.5 and 5 GeV 2 , which is precisely the region 



where most of the recent NMC data on C, Ca and other nuclei have been taken [|T[ 

For light nuclei, such as C, the dominant mechanism for nuclear shadowing involves 
the double scattering of the projectile from two nucleons. Higher order terms (multiple 
rescattering) in the Glauber expansion attenuate the incident flux of vector mesons (or 
of the hadronic state X for the diffractive component |ITJ) as they traverse the nucleus, 
which will be progressively more important as A increases [[|. For the VMD component, 
the magnitude of the attenuation is determined by the mean free path, Ly, of the vector 
meson in the nucleus, Ly = (p uvn)' 1 , where cjvn is the total VN cross section |27j. If 



the mean free path Lx = (p o~xn) 1 of the hadronic state X is independent of the mass 
M x |0, one may take the reabsorption cross section ctxn ~ cryjv ~ 20 — 30 mb, although 



in the Q 2 range covered by the NMC data the F^/F® ratios are not very sensitive to the 
precise value of (Jxn- For the single particle density p in heavy nuclei (A > 16) we use 
the Woods-Saxon (or Fermi) density, while for light nuclei (A < 16) the harmonic oscillator 
(shell model) form is more appropriate [ffijfl . Short range correlations are included through a 
Fermi gas correlation function which puts a "hole" in the two-body density approximately 
0.5 fm wide at 1/2 maximum density ||. The inclusion of correlations has the effect of 
decreasing slightly the amount of shadowing (i.e. increasing the ratio F£/F£) at low Q 2 . 

Having outlined the essential features of the model, we now turn to a detailed comparison 
with the data. In Fig.l we plot the ratio of the structure functions F 2 (normalized to one 
nucleon) for C and D as a function of Q 2 , for various values of x ranging from x = 0.0003 
to x = 0.055. The data represent the complete sample taken by the NMC for C nuclei. 
The overall agreement between the model calculation and the data is clearly excellent. 
In particular, the observed Q 2 dependence of the ratios is certainly compatible with that 
indicated by the NMC data. At large Q 2 (Q 2 > 10 GeV 2 ), the calculated curves are almost 
constant with Q 2 , as would be expected from a partonic (scaling) mechanism fll9|,f29|-|3"2"| . 
However, in the smallest x bins the Q 2 values reach as low as Q 2 ~ 0.05 GeV 2 . Clearly 
this region of Q 2 is inaccessible to any model involving only a partonic mechanism, and it 
is essential here to invoke a non-scaling mechanism such as the VMD model. 



It is common practice in many data analyses to extract slopes in logQ 2 by performing 
simple straight-line fits to the data — such as F^fF® = a + blogQ 2 . To illustrate a 
potential difficulty with using such extracted slopes to discriminate between different models 
of shadowing, we plot in Fig.l (for x = 0.0125) the slope that one would obtain by making 
a linear fit to the logQ 2 variation of F% / F® ■ Evidently the slope thus obtained is negative, 
while the calculated curves have a positive slope. Within the present degree of accuracy, 
the C/D data are equally compatible with a positive or negative slope, and the trend of the 
data as a function of Q 2 clearly favors the positive value obtained in our calculation. 

In Fig. 2 we show the structure function ratio for Ca/D, for x between x = 0.0085 and 
0.07. Again, within the error bars, the calculated x and Q 2 dependence is in excellent agree- 
ment with the data. The apparent small decrease in the ratio for the largest Q 2 points could 
be due to a small nuclear dependence of the ratio of the longitudinal to transverse virtual 
photoabsorption cross sections fl33f , which is assumed to be zero in the data extraction. 
However, the statistics on the data do not allow a definitive statement about this effect. 

In future the NMC will also produce data on the Q 2 dependence of the Sn to C ratio 
P^j j. In Fig. 3 (a) we show the predictions for the structure function ratio as a function of 



Q 2 , for x = 0.0125 (lowest curve), 0.0175, 0.025, 0.035, 0.045 and 0.055 (highest curve). 
Because of the expected higher accuracy of these data linear fits to the logQ 2 de- 



pendence should result in more reliable extracted slopes. In Fig. 3(b) we show the slopes 
b = d (F^/F^ /dlogQ 2 , extracted from the curves in Fig. 3 (a), as a function of x. The solid 
line is the slope for the VMD and iP-exchange mechanisms, while the dashed line represents 
the Pomeron contribution only. We would like to stress that essentially zero slopes would 
be obtained in models where only a partonic mechanism would be utilized |19|p9| -pTl. Ob- 
servation of non-zero slopes would clearly support the hypothesis that the intermediate-Q 2 
region is dominated by non-scaling, higher-twist effects, having a non-trivial Q 2 dependence. 

Data for other nuclei have also been taken recently by the NMC. The A depen- 
dence of the structure function ratio of deuterium, Li, Be, Al, Ca, Fe and Sn to 
carbon has been parametrized as F^ip^/F^ix) oc A a ^ x \ In Fig.4 we show the slope 
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a(x) = d(F 2 A /F?) /dlogA clS cL function of x. The agreement at the low x values is clearly 
very good. At larger values of x (x > 0.07) the data go above unity, and tend to lie slightly 
above the calculated curve. This is a reflection of the fact that the present model has no 
mechanism for antishadowing built in. 

Having obtained excellent fits to the nuclear data, we can be reasonably confident that 
applying the same model to the deuteron ||8| JT^JHB|J5B| ] will yield reliable results. Indeed, 
here one has fewer model parameters to deal with, since all of the shadowing is generated 



through double scattering alone. As was observed in Refs. p,[12|1, the largest uncertainty in 
the calculation of the shadowing in deuterium is the deuteron wave function. Nevertheless, 
the presence of shadowing in the deuteron would be confirmed through observation of a 
deviation from unity in the D/p structure function ratio in the kinematic region where 
Regge theory is expected to be valid. Although the exact value of x below which the proton 
and (free) neutron structure functions become equivalent is not known, it is expected that 
at low enough x, F% — > F£, in which case F^/F^ — > 1 + SF^/F^- In Fig. 5 we show the 
data at very low x taken by the E665 Collaboration M, as well as the earlier NMC data at 



larger x |3? |. The calculated ratio with a small shadowing correction is shown by the solid 
curve, while the result for the case of no shadowing is indicated by the dashed curve. The 
data clearly favor the shadowing scenario. As a fraction of F®, the shadowing correction 
amounts to about 1.5% at x = 10~ 2 up to about 3% for x ^ 10~ 5 . 

In summary, we have seen that the latest data from the NMC and E665 indicate that 
nuclear shadowing in the low- and intermediate-Q 2 regions is controlled by the dynamics of a 
higher- twist, vector meson dominance mechanism. At larger Q 2 this component disappears, 
leaving behind a scaling component which is understood to arise from diffractive scattering 
from the Pomeron component of the nucleon, and which agrees with the approximate Q 2 
independence of the data at large Q 2 . Models based solely on partonic mechanisms can 
therefore provide only limited insight into the physics of nuclear shadowing. 
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FIG. 1. Q 2 variation of the C/D structure function ratio, for various values of x [x = 0.0003, 
0.0055, 0.0125, 0.025, 0.045 and 0.055), compared with the complete NMC data Also shown 

for x = 0.0125 is a linear fit in log Q 2 to the data. 
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FIG. 2. Q 2 variation of the Ca/D structure function ratio, for x = 0.0085, 0.025, 0.035, 0.045, 
0.055 and 0.07, compared with the NMC data [§. 
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FIG. 3. (a) Model prediction for the Q 2 dependence of the structure function ratio of Sn to C, 
for x = 0.0125 (lowest curve), 0.0175, 0.025, 0.035, 0.045 and 0.055 (highest curve), (b) Slope in 
log Q 2 of the Sn/C ratio as a function of x — solid curve is the full result, dashed is the Pomeron 
contribution only. 
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FIG. 4. x dependence of the slope a from the structure function ratio F^/F^ oc A a , compared 
with NMC data on A = D, Li, Be, Al, Ca, Fe and Sn. 
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FIG. 5. x dependence of the D/p structure function ratio, compared with the low-x E665 data 
and NMC data J37| at lar ger x. The dashed curve is the result without any shadowing correction. 
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